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Abstract 

We apply the Color Glass Condensate formalism to photon + hadron production 
cross section in high energy deuteron (proton)-gold collisions at RHIC. We investigate 
the dependence of the production cross section on the angle between the produced 
hadron and photon for various rapidities and transverse momenta. It is shown that 
the angular correlation between the produced hadron and photon is a sensitive probe 
of the saturation dynamics. 



1 Introduction 



The recent RHIC results on the suppression of the hadron transverse momentum spec- 
tra in deuteron-gold collisions, as compared to proton-proton collisions, in the forward 
rapidity region pQ have caused a major excitement in the RHIC community. While the 
observed suppression was qualitatively predicted by the Color Glass Condensate (CGC) 1 
based approaches [3], and verified by more quantitative analysis jU there have been 
newer models introduced recently which can also fit the data [Hj. Measurement of different 
processes, for example prompt photon [7j or heavy quark pair production [S] , at RHIC can 
therefore help establish the dominance of the saturation dynamics in the forward rapidity 
region at RHIC. 

Unlike hadron production, electromagnetic processes such as prompt photon or dilepton 
jH] production do not involve hadronization of the final state and are therefore a cleaner 
process in which to investigate the saturation dynamics. Due to smallness of the elec- 
tromagnetic coupling constant, however, these processes are rare and require high beam 
luminasity and/or long running times of the machine. Nevertheless, in order to clarify the 
underlying physics of particle production in forward rapidity region at RHIC and eventu- 
ally in mid or forward rapidity LHC, it is essential to measure photons, dileptons, etc. at 
RHIC. A measurement of photon production in the forward rapidity region, for example, 
can distinguish between recombination models and the Color Glass Condensate physics 
since one does not expect recombination effects to be relevant for prompt photon produc- 
tion while one expects a suppression pattern for prompt photons very similar to hadrons 
in deuteron (proton)-nucleus collisions in the Color Glass Condensate formalism. 

Two particle correlations are known to be a sensitive probe of saturation dynamics. 
While the two hadron correlation function in mid rapidity is expected to broaden due to 
multiple scattering, it is expected to disappear as one measures two widely separated (in 
rapidity) hadrons due to the small x quantum evolution ^U] . This qualitative expectation 
is verified by preliminary data from RHIC jllj . A more quantitative theoretical analysis of 
the two hadron correlation in the Color Glass Condensate formalism ^2] is however quite 
challenging due to presence of 3 and 4-point functions of Wilson lines and therefore requires 
making approximations (for example, large N c or Gaussian ansatz) and assumptions which 
may affect the outcome of the analysis. Furthermore, since the measured hadrons are at 
low to intermediate transverse momenta, non-perturbative effects may play an important 
role. 

Therefore, in this brief note, we consider production of a hadron and a photon as the 
cleaner and theoretically simpler process in which to investigate angular correlations and 
the role of saturation physics. This process has the further advantage that it is not expected 
to suffer from possible recombination effects and can therefore help establish/constrain 
validity of applying the Color Glass Condensate formalism to deuteron-nucleus collisions 
at RHIC. Experimentally, it is possible to measure this process at RHIC, using the STAR 
detector, for example, even though one will most likely need another deuteron-gold run for 
improved statistics and detector coverage. 

1 See [2] for reviews and extensive references. 
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As shown in one probes the very large x region in the deuteron wave function and 
the very small x region in the target wave function (see Fig. 10 in j3] for the values of x 
contributing to single hadron production). In this kinematics, one can treat the incoming 
deuteron (proton) as a dilute object, consisting of quarks and gluons [T3] with a possible 
nuclear modification of the deuteron wave function which should be very small for the 
typical x and Q 2 involved. On the other hand, the values of x probed in the target wave 
function are very small, 0(1CT 3 — 10~ 4 ). Therefore, it is essential to include the high gluon 
density effects in the target. 

The quark-photon production cross section in quark-nucleus scattering has been calcu- 
lated in the Color Glass Condensate formalism in jS] . We use this cross section and convo- 
lute it with quark (anti-quark) distributions in a deuteron (proton) and quark-hadron frag- 
mentation function to make predictions for hadron-photon angular correlations at RHIC 
and LHC. To probe the smallest x possible in a given process in a hadron-hadron (nucleus- 
nucleus) collider environment, one needs to be in the forward rapidity region. RHIC has 
the further advantage of having a very unique forward rapidity capability, with plans for 
detector upgrades. One therefore can expect to be able to measure the hadron-photon 
correlation in a wide kinematic region at RHIC. For example, the STAR collaboration 
will be able to measure hadrons and photons at rapidities y = and 4 and study their 
correlations. 



2 The Scattering Cross Section 

The scattering cross section for production of a massless on-shell quark with momentum I 
and a real photon with momentum k was derived in jU] . It is given by 

d<7ff (pM^(Q7(k)X e 2 aem k - 1 + (I1)2 



d?b t dk 2 t dl 2 t d Vl dy t d0 ^(2n) 2 k 2 t Js~ [jfe- f t _ l~k t f 

2r k- l t k t cos6 + k~(p- -k~)l 2 + r(p- -r)k; 

dr t r t Jo[r t \h + h\] N(b t , r t , x g ) 



S[x - 4=e M - 

V s V s 



where the incoming quark has momentum p, the photon and outgoing quark rapidities are 
defined via k~ = -^e y "' and l~ = -^e Vl - The angle 6 is the opening angle between the 

final state quark and photon defined as cos 6 = jjf-, with respect to the produced quark 
axis. The dipole cross section iV satisfies the JIMWLK equation and has all the multiple 
scattering and small x evolution effects encoded. It is defined as 

N{b t , r t , Xg ) = ^Tr <l- V\x t )V{y t ) > (2) 

with b t = (xt + yt)/2 and r t = x t — yt- The dipole cross section depends on Bjorken x g 
via the JIMWLK renormalization group equations. In the present case, it is related to the 
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photon and final state quark rapidities and transverse momenta via 



[k t e-y + l t e- 



VI] 



(3) 



In order to compute the hadron + photon production cross section, we would need to 
convolute the above partonic cross section with the quark and anti-quark distributions of a 
deuteron (proton) and quarh-hadron fragmentation function. However, before doing that, 
it is instructive to investigate in some detail, the properties of the above cross section. 
Specifically, we would like to investigate the dependence of this cross section on the angle 
between the produced quark and photon. In order to do this, we isolate the parts of the 
cross section which depend on the angle and ignore the rest of the kinematic factors for 
the moment. Therefore we define the angle dependent part of the cross section as 



1(9) 



[2l~ k~ l t k t cos9 + fc-(p- - k~) l\ + Z~(p- - r) kf 



[k~ h - 1- 



dr t r t J [r t \lt + k t \)N(b t ,r t ,x g ) (4) 



Without making any assumptions about the specific form of the dipole cross section N, 
it is clear from the factor - — =A — ^— that 1(9) diverges when the momenta of produced 

[fe- l t — l~ktY 

quark and photon are parallel [9 — > 0). This is the standard collinear divergence present 
in perturbation theory and is not affected by saturation physics. On the other hand, when 
l t = —kf, then \k + k t \ — > and the integral over the dipole size r t is divergent. This is 
taken care of by demanding color neutrality of the target which would cut the dipole size 
off at r t ~ lfm. It should be noted that this is different from the single hadron production 
case where, at finite transverse momentum of the produced hadron, the integral over dipole 
size r t is always finite. 

To proceed further, we need to know the dipole profile N(b t ,r t ,x g ). It satisfies the 
JIMWLK equation and therefore can be obtained by solving the JIMWLK equation with 
a suitable boundary condition. However, this is highly nontrivial and in practice, it is much 
more convenient to use one of the available parameterizations such as the KKT profile |3] . 
In this model, the dipole profile is given by 



N(b t ,r t ,x g ) 



exp 



\[^Q 2 MV (y ' rt) 



l 



(5) 



where the anomalous dimension j(y,r-t) is 



i(y,r t ) 



and 



1 + 



€(y,r t ) 



tiy,rt) + yj2Z(y,r t ) + 28C(3), 



(6) 



(X/2)(y-y ) 



(7) 



The saturation scale is given by Q s (y) = Qoexp[X(y — yo)/2] with y = lnl/x g and y = 
0.6, A = 0.3. This is the form of the dipole profile we will use for our numerical analysis. 
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However, it is instructive to consider the analytic behavior of 1(0). This is possible in 
two limits; when j(y,rt) is either 1/2 or 1. The former corresponds to the BFKL and 
saturation region while the later is the case when one has the standard DGLAP anomalous 
dimension. It turns out that the case 7 = 1/2 is very close to the actual value of the 
anomalous dimension when one does the integral numerically. Therefore, we use this value 
in our analytic estimates. The integral over the dipole size r t can be then done exactly and 
gives 

- 16 

dr t r t J [r t \l t + k t \]N(b t ,r t ,Xg)\ 1=1 / 2 



16[Zt+fcjPl# 



Strictly speaking, the integral over r t is divergent when = n if \l t \ = \k t \ and needs to 
be regularized by confinement scale r t ~ lfm. To avoid the numerical complications of 
regularization, in this work we stay away from = tt point and keep in mind that the 
expression in (jBJ) is not valid 2 at = tt. Using this expression in 1(0) gives 

1{e) , k t e-y^[l t + k t r 1 1 

1 ) i t [( k -)h-kt] 2 Ql + 1 ) 

This form of 1(0) shows the angular dependence of the cross section most clearly and can 
serve as a guide understanding the exact numerical results shown later. Note that the 
above expression does not depend on whether k t is larger or smaller than the saturation 
scale Q s and is therefore vaild in both the saturation and BFKL regions to the extent that 
7 ~ 1/2 remains a good approximation. 

In Fig. (^) we show the function 1(0) where both the produced photon and quark 
are in mid rapidity (RHIC) and for the case where both have equal transverse momenta. 
The solid line is when both quark and photon have transverse momenta It — kt — lGeV 
and the dashed line is when both transverse momenta are 3GeV. For ease of comparison, 
we normalize the two curves at the smallest angle. At the smaller momentum lGeV, the 
"away" side (0 —* 0) correlation is flat while for higher momentum 3GeV, the away side 
peak is returning. This is mainly due to the increase of x g , as given by (jSJ) and therefore a 
decrease of Q 2 s (x g ) which indicates a weakening of the gluon saturation effects in the target 
nucleus. It should be noted that introduction of an "isolation" procedure experimentally 
in the photon measurement would cutoff the collinear divergence present in 1(0) which is 
apparent in the small angle limit in Fig. Since the isolation cuts needed will depend 
on the detector geometry after a possible upgrade and are not known at the moment, we 
do not consider an isolation cut on the photon. 

In Fig. (J2J) we show the angular correlation function 1(0) at different rapidities. The 
solid line shows the case when both particles are in mid rapidity while the dashed line is 

2 Also, due to the choice of 7 — 1, the high transverse momentum limit of the dipole cross section is not 
reproduced correctly, but this is of minor importance for 9 — > 7r since in this limit the magnitude of the 
argument of the dipole cross section, |Zj + kt\ — ► 0, is always less than the saturation scale. We emphasize 
that in the numerical analysis, the exact form of 7 is kept so that the high transverse momentum limit of 
the dipole cross section is correctly reproduced. 
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Figure 1: The angular part of the cross section, 1(9) in mid rapidity RHIC. 

when the produced quark is at y = 4 while the produced photon is in mid rapidity. In 
both cases, the transverse momenta of both particles is 3GeV. This time, however, the 
two curves are normalized at the highest angle. As before, the angular correlation function 
diverges in mid rapidity when 9 — > due to the unregularized collinear divergence, we 
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Figure 2: Rapidity dependence of 1(9) at different transverse momenta. 

note that there is no collinear divergence (9 — > 0) in the case when the rapidity separation 
is so large, Ay = 4, and when l t = k t = 3GeV while the mid rapidity correlation function 
is collinear divergent. 

Finally, we consider the case when both produced partons are in the forward rapidity 
region, y\ = y 7 = 4. In this kinematics, one is restricted to small momenta due to the 
available phase space so we consider l t = kt = lGeV as well as It = kt = 1.7 'GeV which 
is very close to the upper limit allowed by kinematics. Note that the two curves are not 
normalized like the previous figures. Again in the small angle limit, one is probing the 
collinear divergence present in the angular function which shows up as sharp rise of 1(9). 
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In a realistic experimental set up, this will be tamed by the isolation cut imposed on the 
photon. 




Figure 3: The angular correlation 1(B) in the forward rapidity region y\ = y 1 = 4. 



To relate our results for the parton production to the case hadron + photon produc- 
tion, we need to convolute the partonic cross section in with the parton distribution 
function in a deuteron and the quark-hadron fragmentation function. The hadron +photon 
production cross section can then be written as 



d<J dA^h(q)f(k)X 

d?b t dkt dqt dy 1 dyh d6 



f ^ j dx f q/d (x,Q 2 ) D h/g (z,Q 2 



da q(p)A-,q(l) 1 (k)X 

d?b t dk\ dl 2 dy 1 dyi dd 



(10) 



where a sum over different quark and anti-quark flavors is understood and y^, q t are the 
rapidity and transverse momentum of the produced hadron. We have neglected hadron 
masses so that the rapidity of the produced quark ans hadron are the same. The momentum 
fraction z is defined as z = qt/h — Q~/l~ an d its minimum value in (jlUj) is given by 
z min = ^[ 1 kt ] e yh ~ Vl . Since the nuclear modification (shadowing) of the deuteron wave 

function is expected to be small in this kinematics, it is ignored. We use the GRV98 set 
of parton distribution functions as well the KKP hadron fragmentation functions . 
The integral over x is simple and can be done using the delta function present in the 
partonic cross section so that there is effectively only one integration to perform. 

In Fig. (j3J) we show the invariant cross section given by (fTUj) at mid rapidity RHIC 
for three different transverse momenta; qt = kT = 1.5GeV, qt = k t = 3GeV and qt = 
kt = 5GeV. The lowest transverse momentum of the hadron considered is ~ 1.5GeV 
since the fragmentation functions of KKP start at Qq = 1.5GeV. Again, all the curves 
are normalized to match at the lowest angle for ease of comparison. As is seen, the 
away side correlation is much less for smaller momenta, which is this case is comparable 
to the saturation scale of the nucleus Q s ~ l.bGeV. The largest considered transverse 
momentum, 5GeV, is much larger than the saturation scale and most likely, even larger 
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than the so called extended scaling scale. The cross section is again much larger for the 
collinear configuration of the hadron and photon. This will be modified by the isolation 
cuts imposed by the future detector setup at RHIC. However, the rise of the correlation 
with higher transverse momenta is clear even without any isolation cuts. Furthermore, 
the isolation cuts imposed on collinear photons will not affect the large angle correlation 
even though it will be much easier to judge whether the away side peak is as large as, or 
comparable to, the near side peak which is not possible to tell without an isolation cut. 
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Figure 4: The invariant cross section given by ([10)1 at different transverse momenta. 

In Fig. (jSj) we show the invariant cross section as given by (fTUj) at different rapidities: 
We consider the cases when both hadron and photon are at the same rapidity; either 




Figure 5: The invariant cross section (JTUJl at different rapidities, 
mid rapidity (y = 0) or forward rapidity (y — 4), and the case when they are far apart 
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in mid rapidity, i.e. when the photon is in mid rapidity and the hadron is in the forward 
region. We note that in the kinematics when the two particles have the same transverse 
momenta but are widely separated in rapidity, there is no collinear divergence due to the 
large factor of e Vl ~ Vh in the denominator. Physically, this is easy to understood since 
two particles widely separated in rapidity can not be collinear for comparable transverse 
momenta. When the two produced particles are at the same rapidity, there is always a 
collinear divergence for comparable transverse momenta. We finally note that at rapidity 
of y = 4, we are very close to the edge of the kinematic limit for transverse momenta 
< 2GeV which limits our ability to study the change in the away side peak with increasing 
transverse momenta. Having a more refined hadron or photon detection capabilities in the 
intermediate rapidity region will greatly help with this. 

3 Discussion 

We have calculated the hadron + photon production cross section in high energy deuteron- 
nucleus collisions at RHIC energy and investigated the dependence of the cross section on 
the angle between the produced hadron and photon. It is shown that the magnitude of 
"away" side peak depends sensitively on the transverse momenta of the particles produced 
and can thus be a sensitive measure of the importance of gluon saturation physics. Due 
to the collinear divergence present in the cross section, one needs to introduce an isolation 
cut for the photon. Since this is largely a matter of detector kinematics and capabilities 
which will improve with the possible upgrades at RHIC in the future, we do not include 
isolation cuts and leave this for future studies. 

We used the dipole model of KKT for our quantitative studies. It is probably worthwhile 
to do this analysis using the other parameterizations available, once the experimental setup 
becomes more clear. The KKT dipole profile has been used successfully to fir the forward 
rapidity data without any model assumptions |5]. However, in mid rapidity RHIC, KKT 
dipole profile seem to require a much steeper anomalous dimension JE] or extra model 
assumptions [3]. Therefore, one should be cautious using the KKT parameterization in 
mid rapidity RHIC. 

Furthermore, one can avoid the difficulties of the photon isolation cut procedure, by 
using hadron + dilepton production cross section in order to study the angular correlations 
and the role of saturation dynamics. The dilepton pair invariant mass will regulate the 
divergence of the photon + hadron cross section and therefore is a cleaner process to 
measure even though the rates are smaller. We leave this for future studies. 

Finally, with the expected proton-lead run at LHC coming up in the near future, one 
will have a much broader rapidity coverage which will be helpful in many ways. Due to 
the larger center of mass energy compared to RHIC, the saturation scale will be larger 
at comparable rapidities so that the saturation effects will be stronger, and therefore, 
more robust. Furthermore, some of the proposed detectors such as CMS are expected to 
have very forward measurement capabilities. For example, it may be possible to study 
correlations at LHC between particles separated by 5-7 units of rapidity, and at the same 
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time, have quite a large phase space in transverse momentum available. This will greatly 
help our ability to study the different kinematics regions of the Color Glass Condensate by 
dialing the rapidity and transverse momenta of the produced particles. 
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